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We study a model for anisotropic singlet pairing in .A3C60, using a realistic 
model for the Fermi surface in a hypothesized orientationally-ordered doped crys- 
tal. Anisotropic solutions are studied by combining numerical solutions to the gap 
equation in the low-temperature phase with a Landau expansion near the mean field 
critical temperature. We focus on a class of three-dimensional nodeless d-wave so- 
lutions to the model, which exhibit a fully developed gap everywhere on the Fermi 
surface, but non-BCS temperature dependence of the order parameter in the con- 
densed state, a relatively large value of the low-temperature gap, 2A/kT c , and 
non-BCS structure in the quasiparticle spectrum near the gap. 
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The discovery of superconductivity in the alkali-doped fullerenes [|l| has raised a number 
of interesting questions about the interplay of superconductivity and the unique structural 
properties exhibited by these systems. Exchange of intra- and intermolecular phonons of the 
crystal have both been examined as possible pairing mechanisms for this system |2|-Q . These 
studies have considered a simple pairing scenario in which the condensation is presumed to 
occur on a Fermi surface consisting of a single sheet, with the resulting pair amplitude 
taken to be symmetric under all the point group symmetry operations of the crystal. This 
choice is conceptually simple and furthermore it is apparently motivated by a growing body 
of experimental work on the condensed state which indicates that the gap is nonvanishing 
everywhere on the Fermi surface || . 

In this paper we develop a theory of pairing in A 3 Cqq which treats the structure of the 
Fermi surface of this system more realistically. The model we develop below is actually rele- 
vant to an orientationally-ordered doped system, and in this sense remains highly idealized. 
Nevertheless, we find that even at this level the structure and topology of the three-band 
Fermi surface of this system has important effects on the pairing dynamics. In particular, we 
discuss a class of stable, anisotropic, but nodeless solutions to the gap equation for the three 
band model, in which an "axial" paired state is constructed from a complex combination 
of the anisotropic pairing fields of e g symmetry. Interestingly, this solution does not exhibit 
either the line or point nodes expected of an anisotropic state owing to a very interesting 
interference between this pair field and critical points which occur at intersections between 
the various sheets of the Fermi surface. In addition to breaking the symmetry of the 
high-temperature phase of the doped crystal, this phase is distinguished by a non-BCS tem- 
perature dependence of the order parameter in the condensed state, a relatively large value 
of the low temperature gap, 2A/fcT c , and non-BCS structure in the quasiparticle spectrum 
near the gap. It is an interesting speculation that, unlike the mean field critical temper- 
atures that are predicted for the related fully symmetric a g (s-wave) paired state, these 
results will prove to be relatively resistant to the effects of the short range interelectronic 
repulsive interactions that are omitted from our treatment. However, the fate of this phase 
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in the presence of disorder remains a very open and relevant question. 

Here we briefly present our model Hamiltonian, discuss the symmetry of the low temper- 
ature pair field, discuss our numerical low temperature solutions to the gap equation, and 
develop a Landau theory for the phase transition near the critical temperature. A number 
of details of this treatment will be developed in an expanded discussion to be presented 
elsewhere. 

The conduction band of the alkali-doped fullerenes is derived from an orbitally three 
fold degenerate t\ u state of the isolated molecule. The hopping amplitude between nearest 
neighbor sites can thus be represented by a 3x3 matrix, and to describe it we adopt the 
parameterization of Gelfand and Lu |J , so that H t = ^| [T^z, j)]^j Ua where creates 
an electron with spin a in orbital « on site i. The conduction-band width is then controlled 
by the scale parameter t, while the symmetry of the conduction states is fixed by the form 
of the nearest neighbor 3x3 real matrices T(i,j) given in ||. Energies in this paper will all 
be measured in units of t; the full bandwidth of the conduction band in the orient at ionally 
ordered doped state is W = 42t. 

The Fermi surface (FS) of the model at half filling (three electrons per primitive cell) is 
shown in Fig. 1 [0. The surface consists of two pieces. One surface forms an empty hole 
pocket at the Brillouin-zone center. The second surface is constructed from pairs of parallel 
tubes which run along the Cartesian (001) directions. Symmetry requires that these two 
sheets of FS "touch" at single points along the (111) directions approximately given by 0.52 
(7r/a)(lll); this symmetry will play a central role in the discussion that follows. 

We now expand this model to include an interaction term H int which is decoupled by a 
BCS factorization which yields: 

Hint= E *U^*-k,-a(Ank) + (A^)* „ kflcr D n ^ kl ,_ a - (A^) — (A nU ) (1) 

with (A n k) = ^(^-k/icr-Dn^i/^k^-a-)- Here the creates a Bloch orbital with momentum 
k and the 3x3 matrices D specify the the orbital symmetry of the pair fields that are coupled 
through the interaction potentials U n ; a condensation in the nth channel may occur when 
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U n < 0. For the fully symmetric state with D\ = I, exchange of phonons of A g and H g sym- 
metry both contribute to an attractive potential U\. The electrons can also be paired in an e g 
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doublet which is spanned by two basis functions with _D 2 ,3 = 
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. Basis functions for pairing in a t 2g orbital triplet can also be constructed, but a calculation 
of pairing susceptibility in this channel shows that it is much smaller at half filling and 
we will not consider it further. We will consider only the effect of a static (non-retarded) 
attractive potential, with the one-electron bandwidth then providing a high energy cutoff 
on the theory. 

We begin by studying the second-order pairing susceptibilities d 2 Q/dA a dA^ = {1/\U\ — 
Xa(T)}, calculated in this model as a function of temperature T for the various pair fields 
listed above when the conduction band is half filled. The x's are all logarithmically divergent 
at low temperature. In the on-site approximation, U\ = U 2 = Us, we find that at a given 
temperature Xa is largest for the fully symmetric a g pair field. In the e g channel x is on ly 
slightly smaller so that the form of the paired state for this system will depend sensitively 
on the relative sizes of the of the interaction potentials U n . For a simple model which in- 
cludes only intra-ball scattering and ignores the direct repulsion between electrons, these are 
predicted to be exactly equal for the symmetric and anisotropic channels, so that condensa- 
tion in the symmetric channel should be weakly favored in such a theory. Nonetheless, the 
anisotropic solution is a locally stable solution to the model, and it seems to provide a better 
representation of some aspects of the experimental data. Moreover, it can be dynamically 
stabilized in any theory for which (3/2)|{7i| < \U 2 \, and so we now examine its properties in 
more detail. 

The zero temperature k-integrated quasiparticle spectra near the gap are shown in Fig. 
2. The top panel is obtained for the fully symmetric paired state, and shows a characteristic 



BCS singularity in the spectral density of the form l/^E-zi) 2 - A 2 . The light curve m 
the second panel is obtained for a real anisotropic pairing field of e g symmetry. Here, the 
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quasiparticle spectrum exhibits a pseudogap with the spectral density vanishing as a power 
law at the Fermi energy and peaking at logarithmic critical points located deeper in the band; 
these can already be resolved in these numerical data. However, the "real" gapless state 
considered in Fig. 2(b) is not the lowest-energy paired configuration in the e g manifold; a 
lower free energy is obtained through the axial combination D± = (l/\^2)(D 2 ±iD 3 ). Below 
we show that this combination minimizes the fourth-order terms in a Landau expansion 
of the free energy near the transition temperature, so that the phase transition from the 
normal state occurs directly to the axial phase. This is also established by a numerical 
solution to the gap equations which we have carried out in the low temperature phase, and 
which indicate that the axial phase has a lower free energy than the simple "real" pairing 
states of e g symmetry. The reason for this is apparent in Fig. 2(b) which shows that the axial 
combination (heavy curve) stabilizes a ground state with a well-developed non- vanishing gap 
everywhere on the Fermi surface. 

Nodeless d-wave pairing in a "two-dimensional" system is also possible, and has been 
discussed in the context of high T c superconductivity in the cuprates. For this situation 
the full three-dimensional Fermi surface is presumed to take the topology of a cylinder that 
opens along the k z axis, so that the complex combination A(k) = (k x ± ik y ) 2 yields a non- 
vanishing pair amplitude of (i-wave symmetry which then neatly avoids the point zeroes at 
the poles. The three-dimensional A 3 C 60 system can not make use of this escape route, and 
employs a rather more devious strategy to avoid the point zeroes that would otherwise be 
expected for the axial phase. 

In Fig. 1 we have also displayed the trajectories of the nodes of the order parameter, 
calculated as the zeroes of the diagonal matrix elements (^nk-Da^n-k) for the aih pair state 
on the nth band crossing the Fermi level. On each surface these nodal trajectories are on a 
collision course and are required to meet at each of the high symmetry points where the FS 
crosses the (111) directions. 

Along the (111) directions in k-space, however, the situation is more interesting. Here 
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the one-particle kinetic energy operator H t takes the form H t 



v(k) e(k) v(k) 
v(k) v(k) e(k) 



, so that 
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the t\ u molecular state is broken into a two-fold degenerate A 3 state and a nondegenerate Ai 
state. The crossing of this A 3 doublet through the Fermi surface marks the "contact point" 
between the inner and outer sheets of the Fermi surface. The pairing amplitude at this 
critical point can be written as a real traceless 2x2 matrix for each of the e g components. 
These matrices are not independent since in the A 3 subspace (D 2 ) = (1/2)(I + ia y ){D 3 )(l — 
ic y ). Nonetheless, the eigenvalues of D can only vanish if both the diagonal and off-diagonal 
components go to zero simultaneously exactly at the crossing point. This cannot occur since 
[H t , D] is non-vanishing along the A line, so that H t and D cannot simultaneously be rotated 
to diagonal form. Physically, the off-diagonal components of the pairing tensor dominate at 
the crossing point and provide the required non-vanishing pair amplitude. 

The order parameter of this axial phase exhibits a non-BCS temperature dependence 
below T c . This is shown in Fig. 3(a) where the temperature dependence of the pair amplitude 
(A ± ), and E c , the energy of the logarithmic critical point (which would be readily identified 
in any k-integrated measurement) are plotted as functions of temperature. For comparision, 
the standard weak-coupling BCS prediction is shown as the dashed curve. If we identify the 
logarithmic singularity in these spectra as the "gap" for this system, a relevant parameter 
which can then be extracted from this plot is E c /kT c = 2.3, which is considerably larger 
than the familiar weak coupling BCS prediction of 1.76. In fact particle and hole amplitudes 
along the crucial A line are extremely strongly mixed even out to an energy E\/kT c = 3.08, 
which we believe sets the true energy scale for the pairing in this state. Nevertheless thermal 
excitations should exhibit an activated temperature dependence governed by the minimum 
gap, which is seen to be substantially smaller for this phase. 

Also shown in Fig. 3(a) is the temperature dependence of (A 2 ) and (A 3 ), which have the 
same mean field critical temperature, but are stabilized by larger fourth-order coefficients in 
the Landau expansion, and develop weaker ground-state expectation values at low temper- 
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ature. For comparison in Fig. 3(b) we also display the results of calculations for pairing in 
the fully symmetric channel (this is obtained on the same temperature scale by reducing the 
coupling strength U by 2/3 in the symmetric channel) . Although the transition tempera- 
tures are comparable, (Ai) is found to follow the expected weak-coupling BCS temperature 
dependence very closely. 

The transition from the paramagnetic phase to the axial phase can be studied by de- 
veloping a Landau expansion in the order parameters A2 and A3. To quadratic order 
the free energy of this system takes the form: Q = Qq + (1/|Z7| — x)(|A 2 | 2 + IA3I 2 ) = 
Qq + a[(T — T C )/T C ](|A 2 | 2 + |A 3 | 2 ). The symmetry of the condensed state is then deter- 
mined by the fourth-order invariants in the Landau expansion, which can be grouped as 
VL A = 5f![|A 2 | 2 + |A 3 | 2 ] 2 + 25r 2 | A 2 | 2 | A 3 | 2 + ^ 3 (A|A* 2 + A* 2 A 2 ). We have evaluated each of the 
coefficients g by a perturbation expansion from the high-temperature paramagnetic phase 
carried to fourth order. Stability requires that g% > and cubic symmetry requires that 
92 + 93 — 9i — 0. If we parameterize A 2 = A 3 e^, we see that this potential then takes 
the simple form fi 4 = [| A 2 1 2 + |A 3 | 2 + 2g 3 \ A 2 | 2 | A 3 | 2 (cos 20 — 1), so that the symmetry of 
the condensed state is specified solely by the sign of g 3 . For positive g 3 , we see that Q is 
minimized with |A 2 | = |A 3 | and = ±7r/2, which is exactly the axial phase. Conversely, 
for negative g 3 the system orders with = 0, and in this case the potential is invariant 
under an arbitrary continuous rotation in the A 2 — A 3 plane. (This continuous symmetry is 
not a symmetry of the cubic crystal and is ultimately broken at sixth order in the Landau 
expansion.) Numerically we find that a = 0.017/t, g 1 = 0.0086/t 3 , and 5-3 = 0.0029/t 3 . The 
positive g> 3 indicates that the axial phase is indeed favored over the "gapless" e g phases for 
this system and confirms the results of the low temperature analysis given above. 

It is difficult to make a detailed comparison betwen this simple model and experimental 
data. Nevertheless, we do note that the temperature dependence of the gap has been 
studied by Lieber and coworkers using scanning tunneling spectroscopy and this yields 
a classical mean-field "weak coupling" temperature dependence of the gap in the ordered 
state, but with a non-BCS ratio A/kT c = 2.6. This result has been interpreted as evidence 
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of extremely strong coupling in this system, a result which is surprising (as the authors of 
H note) since this gap is known to be a relatively small fraction of the phonon frequencies 
generally thought to be relevant for scattering the conduction electrons at the FS. Infrared 
reflectance measurements also yield a relatively large gap, A/kT c = 1.5 — 2.5, and this has 
been interpreted as resulting from a distribution of gaps in experimentally inhomogeneous 
disordered materials ||. The temperature dependence of the muon-spin relaxation rate has 
cleanly ruled out an order parameter with nodes on the Fermi surface [||] which, as we have 
seen, is consistent with pairing in either the fully symmetric or axial channels. 

The most provocative property of the axial phase is that it exhibits chiral- symmetry 
breaking in the pair field. This symmetry breaking, if found in this system, should manifest 
itself in any macroscopic property that derives from the low-lying quasiparticle excitations. 
We should point out that for the simple pairing model studied here, the global minimum of 
the free energy is actually found for a condensation in the fully symmetric a g channel. The 
influence of more realistic potentials on this result and of the effect of orientational disorder, 
which seems to be an unavoidable experimental difficulty for the doped fullerenes at present, 
remain to be explored. 

This work was supported by the DOE under Grant 91ER45118, and by the NSF-MRL 
under Grant 91 20668. 
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FIGURES 

FIG. 1. The two sheets of the Fermi Surface in the three-band model of ^Ceo- The outer 
surface is cut away to reveal the inner surface. The lines denote zeroes of the pairing fields for the 
A2 and A3 states (light and heavy curves, respectively). 

FIG. 2. Quasiparticle spectral densities near the gap for pairing in the fully symmetric a g state 
(a), for a real e g state (b, light curve) and for the axial e g state (b, heavy curve). 

FIG. 3. Temperature dependence of the order parameter, and prominent features of the spectral 
densities for pairing states of different symmetries. In (a) the filled circles denote (A±), and the 
filled squares denote E c , the location of the logarithmic peak in the density of states for this 
phase. For comparison,we include the expectation values in the real states (A2) (diamonds) and 
(A3) (open circles). The dashed line is the weak-coupling BCS result normalized to the mean-field 
transition temperature. The chain-dot curves are the results of Landau expansions in the BCS 
model and in our anisotropic model near T c . In (b) the filled circles denote (Ai) for the fully 
symmetric channel, and the filled squares are the locations of the peaks in the spectral density, 
which follow the weak-coupling BCS model (dashed curve) quite closely. 



10 



